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Abstract 

The Bethe-Salpeter (BS) equation in the ladder approximation is 
studied within a scalar theory: two scalar fields (constituents) with 
mass m interacting via an exchange of a scalar field (tieon) with mass 
[i. The BS equation is written in the form of an integral equation in 
the configuration Euclidean x-space with the kernel which for stable 
bound states M < 2m is a self-adjoint positive operator. The solution 
of the BS equation is formulated as a variational problem. The non- 
relativistic limit of the BS equation is considered. The role of so-called 
abnormal states is discussed. 

The analytical form of test functions for which the accuracy of 
calculations of bound state masses is better than 1% (the comparison 
with available numerical calculations is done) is determined. These 
test functions make it possible to calculate analytically vertex func- 
tions describing the interaction of bound states with constituents. 

As a by-product a simple solution of the Wick-Cutkosky model for 
the case of massless bound states is demonstrated. 
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1 Introduction 



The ladder Bethe-Salpeter equation (BS) is an effective if not unique instru- 
ment to study bound states within the framework of quantum field theory 
(see, for example, [1, 2, 3, 4, 5]). However, the exact analytical solution is 
obtained only for the Wick-Cutkosky model in a particular case when the 
masses of intermediate particle and bound state are equal to zero [2, 5]. For 
other cases analytical solutions are not found and numerical methods are 
applied. Numerous literature is devoted to numerical solutions of the ladder 
Bethe-Salpeter equation for different particles, propagators and interactions 
to get acceptable description of spectrum and other characteristics in nuclear 
and particle physics. Although computer numerical methods are very power- 
ful, we are convinced that a simple method to get analytical solutions of the 
Bethe-Salpeter equation in the form of known special functions with accu- 
racy around 0.1-1.0% would be very useful for qualitative and semiquantative 
analysis of physical phenomena, the more if the accuracy can be improved. 
At present, the accuracy around 1% is undoubtedly acceptable in particle 
phenomenology at low energies. 

A generally accepted approach to solve the BS equation ( I do not know 
exceptions ) consists in investigation of this equation in the momentum space 
where the BS-amplitude is expanded on a suitable basis in the Euclidean 
apace R 4 to reduce the original four-dimensional integral equation to an 
infinite set of coupled one-dimensional equations which should be solved by 
numerical methods (see, for example, [2, 6, 7] and many-many other papers). 

The goal of this work is to improve the variational approach to the bound- 
state problem of two scalar particles with mass m interacting through the 
exchange of a scalar particle with mass /i within the ladder Bethe-Salpeter 
equation. In addition, variational calculations are attractive because they 
give restrictions on computable parameters from one definite side, meanwhile 
approximate calculations do not indicate what side we approach to an exact 
value from. Variational methods were applied to the BS equation long time 
ago [8, 9] and now they are used directly in the Hamiltonian formulation of 
quantum field theory (see, for example, [12]). Our idea is the following. In 
the configuration Euclidean x-space the BS equation can be rewritten as an 
integral equation with the kernel which for stable bound states M < 2m is 
a self- adjoint positive symmetric operator. Thus, the problem of solution 
is reduced to searching eigenvalues and eigenfunctions of this kernel. As 
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a result, the Bethe-Salpeter amplitude is expanded over this ortho normal 
system of functions and this series is nothing but the representation of the 
BS-amplitude over the Regge poles in the s-channel. The same representation 
of the BS-amplitude can be written in the t-channel. In this approach, the 
explicit form and normalization of vertex functions is naturally determined 
by the eigenfunctions. Besides, this representation makes clear the role of 
so-called abnormal states. Abnormal states are part of the full orthonormal 
system of functions. They are needed to get the correct representation of 
the Bethe-Salpeter amplitude over bound states. The experimental status of 
abnormal states is not clear up to now. I think this problem deserves special 
investigation. 

Since the kernel is a self-adjoint positive and symmetric operator, the 
task to find eigenvalues and eigenfunctions can be formulated as a variational 
problem. It turned out that quite simple test functions describe the bound 
state masses with accuracy less than 1% (we compare our results with the 
numbers obtained in [6]). 

It is interesting to discuss the connection between the Schrodinger and 
BS equations and to estimate qualitatively the value of the coupling constant 
for which relativistic corrections should be taken into account. It turned out 
that the nonrelativistic limit is realized for very small coupling constants 
a < 1CT 3 . In other words, any situations when binding energy is larger than 
1% of constituent masses, require relativistic description. 

2 The Bethe-Salpeter equation in the self- 
adjoint form and its formal solution 

Within a simple quantum field model we shall consider the simplest case of 
the Yukawa interaction of charged scalar particles (" constituents" ) described 
by the field $ with mass m and neutral intermediate bosons ("tieons") de- 
scribed by the field <fi with mass /i. The Lagrangian density is 

L(x) = $ + (D-m 2 )$ + i0(D-/! 2 )0 + #$ + $0. (1) 

For simplicity, in the following we put m — 1. It means that all dimen- 
sional variables (^-coordinate, energy, momentum and masses) are expressed 
in units of the constituent mass m. This model is frequently used as the 
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simplest pattern of QFT in many discussions, although this system is con- 
sidered to be unstable from a strict point of view because the Hamiltonian 
is not bounded from below (see, for example, [10]). Nevertheless this model 
has been investigated by various methods (see, for example, [2, 1, 11, 12] and 
references therein). 

A two-particle bound state corresponds to a solution of the homogeneous 
Bethe-Salpeter equation for a four-point amplitude. This equation in the 
ladder approximation in the configuration Euclidean x-space looks like 

D(x 1 ,x 2 ; x 3 , x 4 ) = D(xi - x 2 )5(x 1 - x 3 )5(x 2 - x A ) (2) 
+g 2 Jj dz 1 dz 2 D(x 1 - x 2 )S(x 1 - z 1 )S(x 2 - z 2 )D(z 1: z 2 ; x 3 , x 4 ), 

where the propagators of constituents and tieons are 

^^/(2^-l^^^ K ^ < 3 > 

[ ) J (2vr) 4 /i 2 + F (27t)^ l[ ^ h 

Let us write down equation (2) in the center-of-mass frame, i.e. let us intro- 
duce the variables 

y y i y' i y' n y" n y" 

X\ = x+-, x 2 = x--, x 3 = x + — , x A = X - — , z x =x + — , z 2 = x - — . 
Then 

□ (xi, x 2 \ x 3 , x A ) = n(x - x'] y, y') 

and equation (2) reads 

a(x-x';y,y') = y^)6(x - x')S(y - y')^/) (4) 
g 2 JJ dx"dy"D(y)Tl(x - x", y - y")U{x" - x'; y", y'), 



+ 



ii[x — x ,y — y ) — b \ x — x H — \ b I x — x — 
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Now let us go to the Fourier transform 

D(x-x';y,y') = J j^y/^^y') 

and define the function 

n p (y,y') = ^D(i)Z p (y,y')^D{^) 
Then equation (4) becomes 

Z p (y, y') = 5(y - y') + g 2 J dy"K p (y, y")Z p (y", y'), (5) 
with the symmetric kernel 

K p (y,y') = ^D^)U p (y-y')^D(^) (6) 

r dk e ik (y~y') 

r dk e iHy-y') 
" J (2^) 4 (P + l + f) 2 -(A;^' 

We want to emphasize that equation (5) represents the Bethe-Salpeter 
equation in the form of the integral equation with the symmetrical kernel 
K p (y,y'). If we consider the problem of stable bound states, i.e. p 2 = — s = 
—M 2 and M < 2 in (7), then the kernel K p (y,y') is the self-adjoint positive 
symmetrical normal operator on the functional space L 2 . This kernel can be 
represented in the form 

K P (y,y') = J2U Q (y,p)A Q (s)U Q (y,p), (s = -p 2 ) (8) 
Q 

where the functions Aq(s) and Uq{i),p) are the eigenvalues and eigenfunc- 
tions of the kernel K p (y,y'), i.e., 

A Q (s)U Q (y,p) = J dy'K p (y,y')U Q (y',p). (9) 
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Below we discuss the quantum numbers Q. 
The functions {Uq(i/,p)} are orthonormal 

/ dy U Q (y,p)U Q/ (y,p) = 5 QQ ; ^U Q (y,p)U Q (y' ,p) = 5(y - y') 
J Q 

The general properties of the eigenvalues Aq are characterized by the 
traces 




The trace (10) for the kernel (6) is finite T n < oo for n > 3 in the case 
/j, — and for n > 2 in the case /i > 0. It means that eigenvalues decrease 
as Aq(s) ~ with an appropriate constant c. 

If the eigenfunctions and eigenvalues are found, one can write 

5(y - y 1 ) + g 2 K p (y, y 1 ) = ]T U Q (y,p)[l + g 2 A Q (s)]U Q (y',p). 

Q 

Let us represent 

Q 

then equation (9) looks like 

Zq(p) = 1 + g 2 A Q (s)Z Q (p), Z Q (p) = - i— -. 

1 - 9 a q{ s ) 

One gets 

Z P (y,y') = Y, u Q(y>P) - ] _ L , , U Q {y',p). 
q 1 9 A Ql s J 

and 

n p (y,y>) = Y,\fD(y)U Q (y,p) \ U Q (y' ,p)^D(y) (11) 
q 1 9 IV Q\ S ) 
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Finally, the solution of the Bethe-Salpeter equation (2) in the x-space is 
n( Xl ,x 2 ;x 3 ,x 4 ) = a( x - x';y,y') (12) 

The standard Fourier transform of the amplitude □ in the momentum 
space looks like 

J ...J d Xl dx 2 dx 3 dx 4 e ^ix 1+ k 2 x 2 -k 3 x 3 -k 4 x 4 ) n ^ x ^ X2 . x j 
= (2n) 4 5(k 1 + k 2 - k 3 - fc 4 ) □(«,*). 

Here s = —(ki + k 2 ), t = — (ki — k 3 ) 2 , where k±, k 2 and k 3 are the Euclidean 
momenta. Thus, the amplitude □ as the solution of the Bethe-Salpeter 
equation (2) reads 

D(s,t) = j:V Q ( q ,p) 1 ( , V Q (q',p) (13) 

q 1 y iv Qv b ) 

where p = k\ + k 2 , q = k\ — k 2 , q' = k 3 — k 4 and the vertices are 

V Q (q,p) =i l Jdy e-^y/D(y)U Q (y,p) (14) 

The mass of a bound state with quantum numbers Q is defined by the 
equation 

l=g 2 A Q (M 2 Q ), p=(0,iM Q ). (15) 

Formula (13) gives the Regge representation over the poles in the s- 
channel. The same representation is valid in the t-channel. 



3 Variational solution of the Bethe-Salpeter 
equation 

We start with the following observation. The usual variation approach is 
applied to the BS equation written in the well-known standard form (see, for 
example, [8, 9]) 

Lip = XDip, (16) 
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where L and D are self-adjoint positive operators. The result is 



For a test function \l> we get 



A = mm — — — . 17 



(gLg) 

A " Ao " ■ (18) 



However, if we rewrite equation (16) in the form 



\(f) = VDy VD ■ 0, 0=^^, (19) 



the variational problem looks like 



i =mm (*vBif*) =m J*gj^3 (20) 

A (00) V (-0-D-0) V ; 



and for the same test function \1> as in (18) we get 



1 I (tf £>±Dtf ) 

A*AT= <*5») ' A£Al <21) 
One of the corollaries of the Cauchy inequality for a positive operator L is 

(X Y) 2 < (X L X) (YjY) , VX,FG L 2 . 
This inequality for X = ^ and Y = D^J gives 

< ;_ _ and Ai < A (22) 



Thus, the variational method is more preferable to apply to equation (19) 
and use (20). 

Let us come back to equation (9). The kernel K p (y,y') in (6) for stable 
states p 2 = —s = — M 2 , M < 2m is the self-adjoint positive operator. 
The spectrum of this operator is bounded from above. The main preference 
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of the representation (6) is that variational methods can be used to find 
the spectrum and eigenfunctions of this kernel. Therefore, the solution of 
equation (9) can be reduced to the variational problem in the form (20) 

. , , (UqK p Uq) (UqVDYWDUq) 

Aq(s) = ™ x -iww = n c — (« — (23) 

It is convenient to introduce the function 

U Q (y,p) = y/D(y)*Q(y,P)- (24) 

Then formula (23) reads 

(* Q DIl p D* Q ) 
An (s = max — f- — -f- — 25 

where 

(^ Q DU p D^ Q ) = JJdydy' ^ Q (y,p)D(y)U p (y - y')D(y')V Q (y f ,p), 

{y Q DV Q )= Jdy * Q (y,p)D(y)V Q (y,p). 

The vertex (14) is defined by 

V Q (q,p) = i l jdy e-^D(y)^ Q (y,p). (26) 

The representation (25) will be exploited in this paper. 
Another possibility is to introduce the function 

U Q i», P ) = (27) 

-ply) 

Then formula (23) becomes 

Aq{s) - max , x (28) 

This representation is particularly convenient in the case of the Wick-Cutkosky 
model (/i = 0) for which 

D ^ = WrT> and W) = (2n? - y2 - (29) 



The variational problem in the form (25) and (28) is the key object of 
our investigation. 

The mass of a bound state with quantum numbers Q is defined by equa- 
tion (15). The eigenvalue Aq(M) is a monotonically increasing function when 
the mass M increases on the interval M 2 e [0,4m 2 ], i.e. 

A Q (M 2 ) < Aq(M|) if Mi < M 2 . 

For a test function \I/ in (25) we get 

A Q (M 2 ) > Ag ar) (M 2 ). (30) 

It gives 

j 2 = A Q (M 2 ) = Ag ar) ((M^) 2 ) > A^(MJ). 

Thus, the variational calculations in (23,25,28) give the upper estimation for 
bound state masses 

M Q < M { Q ar) . (31) 
If we want to find the coupling constant for a given mass M, then 

1 =A^(M2)<A Q (M2) = 1 

9(var) 9 

and 

9 2 <9lary (32) 
3.1 Classification of states 

Let us consider the classification of states dictated by equation (9). The 
kernel K p (y,y') depends on the direction of the vector p. It is convenient to 
work in the frame p = (0,p^) = (0,iM). 

Let us start with the case M — 0, i.e. p — 0. Then the ladder Bethe- 
Salpeter equation has 0(4) symmetry in the Euclidean configuration space 
and the eigenf unctions look like 

y nKM (y)=T^(y)R nK (y 2 ), (33) 

T ( { ;\(y) = 0(\yn, R nK (y 2 ) = 0(l) for y - 0, 
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where n and k are radial and orbital quantum numbers and T^> Mk (y) are 
spherical harmonics in the space R 4 (see Appendix). All states are degen- 
erated over four-dimensional orbital quantum number k and the number of 
these states is (k + l) 2 for fixed n and k. 

Now let M > 0. Then the Bethe-Salpeter equation possesses 0(3) sym- 
metry and all states are degenerated over three-dimensional orbital quantum 
numbers / which define the spin of corresponding bound states. The degen- 
eration of the state (33) is taken off and the total number of these states can 
be represented 

1=0 

i.e. after removal of degeneration, the state (33) turns into superposition 
of states with three-dimensional orbital momenta I < k. Usually the state 
with I — k is considered as the desired bound state with spin /, but the 
states with / < k, describing the spin /, are called "abnormal" ones and 
are not considered. From the point of wiew of the representation (33) the 
number k characterizes the behaviour of the wave function in the vicinity 
of zero, i.e. ^ n ni{n} ~ \y\ K for y — > 0. Unconditionally the parameter k 
cannot be a normal quantum number, but for numerical calculations with 
a fixed accuracy k is connected with the smallness parameter because in 
real calculations the additional factor \y\ K leads to decreasing of integrals. 
These argumentations are consistent with the standard approach (see, for 
example, [6]), when the BS amplitude is developed over the four-dimensional 
hyperspherical harmonics Z K [ m (x, 0, (ft) and the obtained system of equations 
is truncated over a "nonphysical" quantum number k. 

Thus, we can characterize a bound state by the quantum numbers Q = 
(n, k, /, to) where (/, to) are three-dimensional orbital quantum numbers, n is 
a radial quantum number and k is connected with the behaviour of a wave 
function for \y\ — > 0. The eigenvalues are degenerated over the magnetic 
quantum number to, i.e. 

A Q (s) = A nKl (s). 

On the other hand, one can say that the quantum number k is connected 
with "time" or the fourth coordinate component y\. It has no special name 
and usually these states are called "abnormal". These states disappear in 
the nonrelativistic limit and usually are not considered. However, the rep- 
resentation of the Bethe-Salpeter amplitude (13) is not valid without these 
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states and their role should be clarified. 

We choose the test functions in the form 

V Q (y,p) = V nK i m (y,p) = Y lm (y)W nKl (\y\,yi) (34) 
= Y lm (y)R nKl (\y\,y A )F nKl (\y\,y 4 ). 

Here Y lm (y) = Y lm {ja. y )\y\ l with = y/|y|. Yi m (n y ) = Y lm (Q,4>) is the 
standard three-dimensional spherical function; R nK i(u,v) is a polynomial of 
the degree (n + k) and F nK i(u,v) is a test function. The constant N nK i is 
defined by the normalization condition 



(^ Q D^ QI ) = jdyD(y)^ nKlm (y,p)y 

n'K'l'm' 

The vertex (14) takes the form 

V Q (q,p) = V nKlm (q,p)^i l Jdye-^D(y)Y lm (y)W nKl (\y\,y 4 ) (36) 
= Yi m (q)V nK i(q,p). 

where 

V nKl (q,p) = JJ dydy 4 {qMP^n^e'^^D (^y 2 + y\) W nKl {\y\^) 
Finally, the BS amplitude (13) looks like 

D(MH S^ v; " ,(, ' p) r^bM K " ,(9>) (37) 

The test function F(u, v) contains variational parameters. We choose the 
test function in the following simplest form 



W(y) = W(\y\, \y 4 \) = e^V^I, (38) 

dk e iky 



D 3 (y) = D 3 (\y\,\y 4 \) 



(2tt) 4 (a 2 + k 2 + bkj) 



2\3 ' 



where a and b are variational parameters. 

Polynomials R nK i(u,v) have to provide the orthogonality condition (35). 
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4 Nonrelativistic limit 



4.1 The "Bethe-Salpeter form" of the nonrelativistic 
Schrodinger equation 

First of all let us represent the nonrelativistic Schrodinger equation in the 
"Bethe-Salpeter form". Bound states of two particles with masses m — 1 
interacting via an attractive potential are described by the equation 

Hip(x) = (p 2 - aV(r)) ip(x) = -e^(x). (39) 

where a is the coupling constant, e is the binding energy and 

x G R 3 , r = V^, p 2 = -<9 X 2 . 

The function V(r) > is supposed to be positive. For the Yukawa potential 
V(r) = ^f- 

Let us rewrite equation (39) in the form 

(p 2 + e) V>(x) = aV(r)^(x), ^(x) = a— V(r) ■ ipU) 

v ' p + e 

and let us introduce the function 

0(x) = sJV(r)^{yL) 

Then we get 

0(x) = a^V{r) ■ ■ ^(x) =ajdx! 5 e (x, x')0(x'). 

Here the kernel S^x, x') for e > is the self-adjoint symmetric positive 
operator 

S e (x,x') = V^.--!— 5(x-x')- v^ 7 ) 

c x -+" e 

= v/y(0 • n e (x - x') • Jvyj, 

= /(2^^rr = s'T^- (40) 
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The solution of the Schrodinger equation (39) is equivalent to the variational 
problem which can be written in different forms 



1 = aA(e) = a max . , - = a max — ; — - — = a max r-. 

<t> (00) V tyVtl)) * 

(41) 

We will show that in nonrelativistic limit M — » 2m the variational problem 
(25) is reduced to (41). 

4.2 The nonrelativistic limit in the Bethe-Salpeter equa- 
tion 

Now let us come back to the Bethe-Salpeter equation in the form (28) and 
consider the polarization operator H p (y — y') 

r dk r dk 4 e -*(y-y')-*k4(w-»i) 

where p = (0, iM) and (pk) = iMk 4 . 
For a state ty(y,p) we have and 



Here A = 1 - and 



dk rdk d ^ 2 (k,k 



(27t) 3 7 2tt (k 2 + k 2 + A) 2 + M 2 ki 



f(k,A; 4 ) = jdy e^D{y)^{y,y 4 ). 
In the nonrelativistic limit 



M 2 

M ^ 2, A = 1 - — -> 0. (42) 



Introducing new variables k = v A q and k 4 = Aq 4 one can get in the limit 
(42) 



P ' J (2tt)3 7 2vr (q 2 + Ag 2 + l)M 



+ 4g 4 2 
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/a f dq - f dq4 ^(v^q^Q) 

~* J (2tt) 3 J 2tt (q 2 + l) 2 + 4g 2 

1 r dq, * 2 (-s/Aq,0) 1 r dk ^ 2 (k, 
4 7 (2tt) 3 q 2 + l ~ 4 y (2tt) 3 k 2 + 



0) 
A 



The test function is chosen in the form ty(y,p) = ^(y, 0) = ip(y). It leads to 



R(p,0) = j d yjdy 4 e-^D(y)^(y) = j-jdy- 



-»\y\ 



-»py 



y 



^(y) 



— / dye-^V(\ymy) 



and 



(*DV) = jj dydy A D(y)^ 2 (y) = —J dy^^ 2 (y) = — (ipVip) 
The binding energy equals 



(2 — M) • ( 1 — 



M 



4 



Collecting all these formulas we obtain (41) 



1 = A(e) = a max — — 
where a is connected with the coupling constant g in the following way: 



a 



9 Z 



9 Z 



167T IQirm 2 



A 



A 



47rm 2 



(43) 



The main result for the solution of the variational problem of the rela- 
tivistic Bethe-Salpeter equation (25) for small binding energy e 1 is that 
it is convenient to choose the test function ^ in the form 



(44) 



*(p) = ^(VAap, Abp A ) 
where a and b are variational parameters. 



15 



5 The Wick-Cutkosky model fi = 

5.1 The mass Mq = in the Wick-Cutkosky model 

As a by-product of the representation (8) equation (9) can be easily solved 
in the case Mq = (see [1, 5]). Equation (9) reads 

U = g 2 AU = g 2 KU 



or 



- ih^l^-^^f^ (45) 

We look for the solution in the form U(y) = \fy I Equation (45) reads 

y 2 $(y) = j£y 2 J dy'U (y - y>)<S>(y>) 
and in the momentum space 

- D »*<*> = T-(i?TTF* ( *>' a = W« < 46 » 

The states are degenerated over quantum numbers m and k — I. Let us 
extract the angular variables 

Hk)=T£... tll (k)A(k 2 ) 

where (k) is the spherical function in the four-dimension space (see 

Appendix) and satisfies the correlations (60) 

^;,,w=o. (*^)^«(*)=^?... w (*)- 

Then equation (46) reads with k 2 = u G [0, oo] 

u A"(u) + (2 + I) A\u) + - -. A(u) = 0. 

7T [U + l) z 
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Let us introduce new variables 



t = 



u-l 
u + 1 



G[-l,l], A(n) = (l-t) l+1 B(t). 



We get 

(1 - t 2 ) B"(t) - 2(2 + I) tB'{t) + - I 2 - 3/ - 2) B(t) = (47) 

The solution is the Jacobi polynomial B(t) = P( l+1,l+1 \t) with 

- - I 2 - 31 - 2 = n(n + 21 + 3) 



5.2 Variation calculations 

For the Wick-Cutcosky model we perform variational calculations for the 
lowest and first excited states and compare our numerical results with the 
numbers obtained in [6]. 

5.2.1 The lowest state 

For the lowest state the test function is chosen in the form 



7T 



and finally 



n(n + I + l)(n + 1 + 2) 



(48) 



MV,P) = W(\y\, M) = e-V^f (49) 



where a and b are variational parameters. Then 



V (k,p) = J dy D{y)W{\yl\y A \)e^ k y = J 




ab + \/a 2 b + ub + v 



[u + v + a?{l + b) + 2aVa 2 b + ub + v] Va 2 b + ub + v 
= R(u,v;a,b), u = k 2 , v — k\\ (50) 



and 
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The eigenvalue A(M) is denned by 



> , <■ ' / „ ,2?, 4o 2 (l + Vb)R 2 (k 2 , v 2 ;a,b) . . 

3 2 A(M 2 ) = « max— \ dk k 2 \ dv n 9 ^ - ^ , v ^ ^,, (51) 



o 



o 



(P + t; 2 + A) 2 + 4(l- A)w 2V 



The value of M is defined by equation (15). It is instructive to compare the 
Bethe-Salpeter mass (51) with the bound state mass of the nonrelativistic 
Schrodinger equation for the Coulomb potential 



e = 2 - M 



a' 
T' 



1 

a 



1 



OL, 



= 1.33... • 10~ 5 , a em = 



8 (l - f ) asch 
1 



137 



The numerical results are represented in Table 1 where we have used a = 
v^Aa-i and b = Abi. 

Table 1. The lowest state in the Wick-Cutkosky model. The coupling 
constant A = 4a for different masses Mj m and comparison 1/a with the 

nonrelativistic value 1/asch- 



M 
m 


e 
2m 


1 


i 

a 


ai 


bi 


A = 4a 


At, [6] 


0.000 


1.00 


0.3535 


0.1586 


0.76 


1.0 


25.21 


25.13 


0.500 


0.75 


0.4082 


0.1671 


0.78 


1.0 


23.93 




1.000 


0.50 


0.5000 


0.1994 


0.80 


1.0 


20.06 


20.01 


1.500 


0.25 


0.7071 


0.3023 


0.84 


1.3 


13.23 




1.900 


0.05 


1.5811 


0.8901 


0.88 


2.1 


4.493 


4.483 


1.999 


0.0005 


15.811 


13.993 


1.00 


4.5 


0.285 


0.285 


O a em 
Z 4 


a 2 
8 


137 


133.9 


1.0 


6.0 


0.0299 





One can see that the variational calculations practically coincide with the 
numerical results obtained in [6]. 

It is important to remark that the nonrelativistic regime, when the Schrodinger 
equation can be used, is realized for very small binding energies e < 10~ 4 . 



18 



5.2.2 The first orbital exited states 

The test functions for the "normal" state with k — 0, 1 = 1 and for "abnor- 
mal" state with k—1, 1 = look like 

& n \y,p) = y W(\y\, \y 4 \), ¥ a \y,p) = y 4 W(\y\, |y 4 |) (52) 

where the test function W is the first function in (38). One can get after 
simple calculations 



V (B) (fc,p) = iJdyyD(y)W(\y\,\y 4 \)e- ik « 



d 

— 2k— R(u, v; a, b) = —2kR u {u, v; a, 6); 



VM(k,p) = Jdyy 4 D(y)W(\y\,\y 4 \)e- 



iky 



d ~ 

= — 2k 4 — R(u, v: a, b) = —2k 4 R v (u^ v: a, b), 
dv 

where the function R is defined by (50), and 

1 ] J (2tt)V 16a 4 (l + v/6) 2 ' 

-- 3 1 



1 J ~" 7 (2vr)V 



(2tt)V 16a 4 ^(l + V^) 2 ' 

The eigenvalue A(M 2 ) is defined by 

2,(n )/nr , 8 7, 1l2 f, 16a 4 (l + Vb)Hk 2 R 2 u (k 2 ,v 2 ;a,b) 
g 2 A (n) (M) = a max— dk k 2 dv — — 2 , 

a ' b 71 { { 3{2 + Vb) (k 2 + v 2 + 1 - + M 2 v 2 

g 2 A (a) (M) = a max— dk k 2 dv - - — - — ^ — - — (53) 

a ' 6 { { 3(k 2 + v 2 + l-^) + M 2 v 2 

The value of M is defined by equation (15). The nonrelativistic binding 
energy of the bound state with / = 1 is 

® 2 i 1 1 1 1 
e = 2 — M = — , and - = = = , a em = — 

16 a / 32 A _ m\ a5 cfe 137 
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The numerical results are represented in Table 2 where a = \fKa\ and b = 
Ab±. Here we should remark that the calculations in [6] contain a mistake 
for I > 1. For the state with / = 1 factor \ should be introduced to get the 
correct value. 

Table 2. The first excited states in the Wick-Cutkosky model. The coupling 
constant A = 4a for different masses Mj m and comparison 1/a with the 

nonrelativistic value 1/agch- 





Vector particles I = 1 


"Abnormal state" 


M 
m 


e 
2m 


1 

&Sch 


a 






= 4a 


X T j 


A, 


ai 


bi 


0.000 
0.500 
1.000 
1.500 
1.900 
1.999 

Z 16 


1.00 
0.75 
0.50 
0.25 
0.05 
0.0005 

a 2 
32 


0.176 
0.204 
0.250 
0.353 
0.790 
7.906 
137 


.0529 
0.0559 
0.0676 
0.1059 
0.3409 
6.4894 

133.9 


0.80 
0.85 
0.81 
0.87 
0.89 
0.95 
0.986 


1.00 
1.00 
1.20 
1.48 
2.56 
5.00 
10.0 


75.62 
71.47 
59.19 
37.77 
11.73 
0.616 
0.0298 


80 
63.5 

16 

5 


75.55 
72.78 
64.27 
48.55 
25.63 
8.98 


0.85 
0.85 
0.90 
0.94 
1.05 
1.30 


1.0 
1.0 
1.1 
1.4 
2.85 
14.0 



One can see that "abnormal states" have practically the same mass as 
"normal" states in the strong coupling regime when the mass defect is large. 
If the coupling constant is small the "abnormal" masses are much more then 
"normal" masses. In addition we can repeat that the nonrelativistic regime 
is realized for very small binding energies. 



6 Massive tieons \i > 

In this section we perform variational calculations in the case fi > for the 
lowest and the first excited states and compare our numerical results with 
the numbers obtained in [6]. 
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6.1 The lowest state 

For the lowest state we choose the test function in the form 

dk * iky 



/dk (D l ^y 
' (^TWTbW (54) 



Here a and b are variational parameters. 
The direct calculations give 



= jdy D(y)D 2 3 (y, a, b) = — ■ a, b), (55) 
1 

Z(fi,a,b) = J J 



dtdu uH(l — t) 



b(bu + 1 - u) {a 2 u + /i 2 (l - u)t(l - t)) 3 



J dy D(y)^(y)e lk y = J dyD(y)D 3 (y, a)e iky = — • V(r, 6, a, b), (56) 



2(4vr 

V(r,9,a,b) 1 <U11 1 



Vbt + l-t [ fl2f + ^ _ t ) + t(1 _ t)r 2 ( cos 2(^) + gg.) 



2 ' 



where we put |k| 2 = r 2 cos 2 (6>) and k\ = r 2 sin 2 (#). 

The mass of the lowest state is defined by the variational equation (25). 
Our formulation of the problem is to find the coupling constant A = 4a for 
given u = — and M — —, then we get 

& " m TO ' & 

. n 2 Z(fi,a,b) 

* = mm— — — — (57) 

a,b J(ji,M,a,b) 

where 

7T 

r 3 cos 2 {6) V 2 {r,6,a,b) 



J(H, M, a,b) = I dr j d6- 



o o 



(r 2 + l-^f) 2 + M 2 r 2 sin 2 (^) 



The numerical results are listed in Table3. 

Table 3. The lowest states in the case \x > 0. 
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m 


M 
m 


\ = 4a 


a 

6=1 


X = 4a 


& 


b 


Xxj 


1.000 


0.000 


42.960 


1.27 


42.960 


1.27 


1.00 


42.96 


1 000 

X .www 


1 000 

x .www 


36 945 


1 1 5 

X . X <u 


36 940 


1 27 


1 10 

X . X w 


36 94 


1 000 

X .www 


1 Q00 

X . <_/ WW 


1 7 293 


66 

w . w w 


1 7 252 


70 

W . 1 w 


1 87 

X . (J 1 


1 7 23 


1 000 

X .www 




1 294 


41 

W .t: X 


1 273 


44 

W .tit: 


2 35 


1 25 

X W . t-i\J 


0.100 


0.000 


25.802 


1.02 


25.802 


1.02 


1.00 


25.80 


0.100 


1.000 


20.692 


0.89 


20.680 


0.91 


1.20 


20.68 


0.100 


1.900 


5.3843 


3.50 


5.262 


0.38 


2.22 


5.227 


0.100 


1.999 


1.0789 


0.089 


1.054 


0.089 


9.00 


1.043 


0.001 


0.000 


25.133 


1.00 


25.133 


1.00 


1.00 


25.13 


0.001 


1.000 


20.023 


0.87 


20.013 


0.87 


1.20 


20.01 


0.001 


1.900 


4.6904 


0.32 


4.501 


0.35 


4.20 


4.483 


0.001 


1.999 


0.3523 


0.031 


0.2900 


0.04 


150 


0.289 



In Table 4 we listed the results of calculations for two cases: (1) one 
parameter a with 6=1 and (2) two parameters a and b. One can see that 
only for small binding energies there is a remarkable difference in numbers. 
Our numbers are quite close to the numerical results obtained in [6]. 



6.2 The first excited states 

For the first excited states we choose the test functions for "normal" (n) and 
"abnormal" (a) states in the form 

*^>-^-«-{te: y ^) (58) 

Repeating as before all calculations we reduce the problem to find the 
coupling constant X = Aa for given u = and M = — to the variational 
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task 



Here 



. TT 2 Zj(fi,a,b) 

A = mm - . — - — , (i—n,a). 

a,b Jj(fi,M,a,by y ' 



Z n (fi, a, b) = 



Z a (n, a,b) = JJ 



dtdu 18w 4 (l -u)t 2 (l -tf 



b(bu + 1 - u){a 2 u + /i 2 (l - u)t{\ - t)) 4 
dtdu 6w 4 (l - u)t 2 {l - t) 2 



{ y/b(bu + 1 - u) 3 {a 2 u + /i 2 (l - u)t(l - t)Y ' 



(59) 



Jj (/i, M, a, b) 



r 3 cos 2 (9)V 2 (r,6,a,b) 



dr / d9- 

iJ (r 2 + l-^f) + M 2 r 2 sin 2 (9) 



l 



K(r,^,a,6) 



4r cos(^) 



+ 1 - t h(/j,, r, 9, t, a, b) ' 



/■ dt t 3 Vb Arsm(U 

V a (r, 9, a,b) = / . = • — — -, 

{ J(bt + 1 - t 3 h(fi,r,9,t,a,b) 



4r sin(^) 



r, 6 1 , t, a, 6) 



a 2 t + /x 2 (l - t) + t(l - t)r 2 cos 2 (9) + 



sin 



bt + l-t / 

In Table 4 the results of calculations for "normal" states with spin I = 1 
and for " abnormal" states with spin I = are listed. The column \ T j contains 
the results obtained in [6] with factor \ to get the correct numbers. One can 
see that our numbers are quite close to them. It should be noted that in the 
strong coupling regime the "abnormal" states begin to play an essential role 
in the representation of the BS amplitude (13). 



Acknowledgments. I would like to thank S.Dorkin for useful discussions 
and help in numerical calculations. This work was supported by the Russian 
Foundation for Basic Research under grant No. 01-02-17200-a. 
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Table 4. The first excited states in the case ji > 0. 





Vector particles / 


= 1 


" Abnormal state" 


rn 


M 
m 


\ n = 4a 


(X 


b 




\ = 4a 


a 


b 


1.000 


0.000 


205.95 


1.6 


1.0 




205.95 


1.6 


1.0 


1.000 


1.000 


183.16 


1.5 


1.1 




195.35 


1.5 


1.1 


1.000 


1.900 


110.86 


1.1 


1.5 




166.26 


1.3 


1.3 


1.000 


1.999 


93.48 


0.9 


1.9 




161.66 


1.3 


1.5 


0.100 


0.000 


80.35 


1.1 


1.1 


80.22 


80.30 


1.0 


1.0 


0.100 


1.000 


63.77 


0.9 


1.2 


63.77 


63.74 


0.9 


1.1 


0.100 


1.900 


16.05 


0.4 


3.2 


15.99 


34.78 


0.6 


3.2 


0.100 


1.999 


4.87 


0.2 


6.2 


4.77 


27.44 


0.4 


4.4 


0.001 


0.000 


75.40 


1.0 


1.1 




75.40 


1.0 


1.0 


0.001 


1.000 


59.06 


0.9 


1.2 




64.24 


0.9 


1.2 


0.001 


1.900 


11.84 


0.3 


3.2 




25.75 


0.4 


4.1 


0.001 


1.999 


1.03 


0.1 


5.2 




10.11 


0.05 


15.0 



7 Appendix. Angular polynomials T^(k) 

The angular polynomials T^(k) = T^ ^k) are symmetric for /ij # 
and 

T w (k) = 

n,n,H3,...,m\ J 

The recurrence relation is 

rfi(*) = }p(l|2...0fciT£ 1) (fc) - ^f^PW.J)5 l2 Tt?{k) 
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In particular 

= 1, T{ 1] (k) = k u 

Tg\k) = k 1 k 2 -^5 12 k 2 , 
k 2 

Tnlik) = hk 2 k 3 - — (M23 + k 2 5 31 + k 3 5 12 ), 
6 

The normalization condition is 

where C}(t) is the Gegenbauer polynomial. 
These angular functions satisfy the formulas 

^l, M {k) = 0, [k^j TSl.^k) = I T^(k). (60) 
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